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The Andreev reflection ampfltude at a clean interface between a half-metaUic ferromagnet (H) 
and a superconductor (S) for which the half metal's magnetization has a gradient perpendicular to 
the interface is proportional to the excitation energy e and vanishes at e = [Beri et ai, Phys. Rev. 
B 79, 024517 (2009)]. Here we show that the presence of impurities at or in the immediate vicinity 
of the HS interface leads to a finite Andreev reflection amplitude at e = 0. This impurity-assisted 
Andreev reflection dominates the low-bias conductance of a HS junction and the Josephson current 
of an SHS junction in the long-junction limit. 
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I. INTRODUCTION 

The experimental observation of a sizable super- 
current through a Josephson junction containing the 
half-metallic ferromagnet Cr02^-''^ has renewed theo- 
retical interest in the superconductor proximity effect 
in half metals-l^^^ Because of t he ab sence of minor- 
ity spin carriers in a half metaljI^^lSl the induced su- 
perconducting correlations must be of the spin-triplet 
type, even if they arise from proximity to a spin-singlet 
superconductor .I^StiSI 

The spin-triplet proximity effect is mediated by An- 
dreev reflectionj25 of majority electrons into majority 
holes. Since such Andreev reflections violate spin conser- 
vation — spin-conserving Andreev reflection reflects ma- 
jority electrons into minority holes — , they require break- 
ing of the spin-rotation symmetry around the half metal's 
spin quantization axis. A prominent mechanism for the 
breaking of the spin-rotation symmetry is a magnetiza- 
tion gradient perpendicular to the interface between the 
superconductor and the half metal, see Fig. [T] which oc- 
curs naturally in the presence of a magnetic interface 
anisotropy that favors a different magnetization direc- 
tion than the magnet's bulk anisotropy.^^ In the litera- 
ture, an interface with a different magnetization direction 
than the bulk magnetization direction is referred to as 
"spin active" W Artificially created spin active interfaces 
have been shown to be responsible for the observation 
of the spin-triplet proximity effect in standard (non half- 
metallic) ferromagnets '^"^ In view of the observe d bi- 
axial magnetic anisotropy of Cr02 thin films ,'^^^^'261 spln- 
active interface have also been proposed as an explana- 
tion for the observed Josephson effect in the Cr02-based 
Josephson junctions,'^ but there is no direct experimental 
evidence confirming this link.'^ 

What makes the spin-triplet proximity effect in a half 
metal particularly interesting from a theoretical point of 
view, is that there is a symmetry argument that strongly 
limits the effectiveness of a spin-active interface as a 
source of spin-triplet superconducting correlations.'^'^ 
This symmetry argument applies to the case of a half 



metal, but not to a standard ferromagnet (with incom- 
plete spin polarization). It essentially forbids Andreev 
reflection of majority electrons into majority holes or vice 
versa if five key conditions are met.'^ Particle-hole sym- 
metry, translation symmetry along the superconductor 
interface, 7r-rotation symmetry around an axis normal to 
the superconductor interface, quasiparticle conservation, 
and the absence of minority carriers. All five conditions 
are met for Andreev reflection of carriers at the Fermi 
level in a clean spin-active interface between a half metal 
(H) and a superconductor (S). Away from the Fermi level, 
particle-hole degeneracy is lifted, and generically the An- 
dreev reflectio n ampli tudes are proportional to the exci- 
tation energy e.EESMlThis energy dependence leads to an 
interface conductance G oc (eV)^, where V is the apphed 
bias, and to a Josephson current / oc max.{kBT, E^)^ 
for a clean SHS junction at temperature T in the "long 
junction" limit (gap Aq mu ch larger than the junction's 
Thouless energy ig^) j '^' l io izH Jn contrast, for junctions in- 
volving a standard ferromagnet, the zero-b ias co nduc- 
tance G is finite, whereas / oc max{kBT, E^)^^^^^ 

In order to open up the possibility of Andreev reflection 
at the Fermi level, one of the remaining four conditions 
listed above has to be lifted. Several options have been 
investigated in the literature, and have been shown to 
lead to a flnite zero-bias conductance G of an HS junc- 
tion and to a temperature dependence of the Joseph- 
son current of a ballistic SHS junction that closely re- 
sembles junctions with a standard ferromagnet instead 
of a half-metal. These include the breaking of the ro- 
tation symmetry around an axis normal to the super- 
conductor interface by a magnetization gradient paral- 
lel to the interface® ^'^ or by spin-orbit coupling in the 
superconductor, or the inelastic scattering of quasipar- 
ticles in the half metal, which effectively lifts quasiparti- 
cle conservation.'^^ 

In this article, we investigate the lifting of the trans- 
lation symmetry along the interface by scattering from 
non-magnetic impurities. This is particularly relevant 
for Cr02-based junctions, as interfaces between the 
metastable compound Cr02 and other materials are no- 



2 



z=0 



t t t t t t 
^ ^ ^ ^ ^ ^ 
^ ^ ^ ^ ^ ^ 



H 



FIG. 1: HS interface with a spin-active interface, originating 
from a magnetization gradient perpendicular to the interface. 
Arrows indicate the magnetization direction. 



toriously poorly defined.'^ We confirm that impurity 
scattering, too, is a viable mechanism for Andreev re- 
flection of majority electrons into majority holes at a HS 
junction at the Fermi level. Moreover, we show that this 
impurity- assisted Andreev reflection remains coherent, 
with a well-defined magnitude and phase, after perform- 
ing an average over different impurity configurations. As 
a result, impurity-assisted Andreev reflection not only 
gives rise to a finite zero-bias conductance, but also to 
a strong enhancement of the Josephson current at low 
temperatures. Our analytical calculations are consistent 
with a numerical analysis of the triplet proximity effe ct in 
a strongly disordered half metal by Asano el al. who 
did not report any suppression of the induced supercon- 
ducting correlations for energies near the Fermi level. 

The detailed outline of this article is as follows: In sec- 
tion [llj we describe the model Hamiltonian used in our 
calculations and review the basic symmetry relations for 
the scattering matrix of an HS interface. Then, in Sec. Ill 
we calculate the Andreev amplitude for an HS interface 
with a single impurity, up to quadratic order in the impu- 
rity potential. We apply our results to an interface with 
a finite density of impurities in Sec. |IV( where we show 
that the presence of impurities at the interface leads to a 
finite interface conductance at zero bias and to a signifi- 
cantly enhanced Josephson current at low temperatures. 
We conclude in Sec. |V| The appendices contain various 
additional results for scenarios not covered in the main 
text. 



II. HALF-METAL-SUPERCONDUCTOR 
INTERFACE WITH PERPENDICULAR 
MAGNETIZATION GRADIENT 

A. Model Hamiltonian without perturbations 

Our calculation builds on the calculation of the 
Andreev reflection amplitudes for a clean half-metal- 



superconductor interface with a perpendicular magne- 
tization gradient by Kupferschmidt and one of the 
authors.!^ Following Ref.[Tni we choose coordinates, such 
that the superconductor occupies the half space z > 
and the HS interface is in the xy plane, see Fig. [l] Pe- 
riodic boundary conditions are applied in the x and y 
directions, with periods and Wj,, respectively. Quasi- 
particle excitations near the HS interface are described 
by the Bogoliubov-de Gennes equatiorPSl 



H^-W e^-W, H = Ho + V, 



where 



Ho 



Ho ^A(r)^72 



and the four-component spinor 



(1) 



(2) 



(3) 



consists of wavefunctions Ucr(r) for the electron and Wo-(r) 
for the hole degrees of freedom. The 4x4 matrix oper- 
ator Ho describes the HS junction in the absence of a 
magnetization gradient and impurity scattering. These 
two effects are described by the perturbation V and will 
be discussed in the next subsection. 

The superconducting order parameter 
A(r) = Aoe*'^e(z), where 0(z) = 1 ff z > and 
otherwise. This step function model is a good 
approximation for tunnelling interfaces of s-wave 
superconductors.!^ For the single-particle Hamiltonian 
Ho we take 



fi2 



'^jJa{z)Pcr + tiw5 (z) , 



(4) 



where m is the effective electron mass (taken to be equal 
on both sides of the interface), and 



^HfT if 2 < 0, 
if z > 0, 



(5) 



with a and the potentials /iH-f, A'h^, and /is repre- 

senting the combined effect of the chemical potential and 
band offsets for the majority and minority electrons in 
the half metal and for the superconductor, respectively, 
and where w sets the strength of a delta-function poten- 
tial barrier at the interface. The operators 

A = ^ + • 6", h = \- • ^, (6) 

project onto the majority and minority components, re- 
spectively, where we have taken 63 to be the unit vector 
pointing along the magnetization direction in the half 
metal. 

The potentials /^Ht: Mh^, and /^s are such that /inti 
/is > 0, and /iH^ < 0. As a result, majority states in the 
half metal with uniform m(r) and states in the normal 
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state of the superconductor are propagating states, with 
Fermi wavenumbers 

fct = ^v/2m/iHt, ks = ^v^2rn^g, (7) 

respectively. The corresponding Fermi velocities are 
= hk^/m and = hko^/m, respectively. Minority 
states in the half metal are evanescent with wavefunc- 
tion decay rate 



1 



2m|/iH;|- 



(8) 



The Andreev approximation Aq <C /is is used throughout 
our calculation. 

The propagation of electrons with Hamiltonian ([T]) is 
described by the 4x4 matrix Green function tj(e;r,r'), 
which is a solution of the Gorkov equation 

(e-H)g(e;r,r') =<5(r-r'). (9) 
The Green function can be written 

g{e- r, r') = ^-(^5 ' ^' ^Oe'''"-^'"""'""^ (10) 

where ky = k^e^ + kyBy is a wavevector parallel to the 
superconductor interface and ry = xe^ + yey.^^ Solution 
of Eq. ^ then gives the result 

g,(£;k||;z,z') = l(t/^,t/^,Vt,n)> (11) 

where the first two column vectors and have the 
form 



1 

1 



/ ^ik^z{e)\z-z'\ _|_ ^^j-^-jg-ifct^(e)(z+z') \ 








(12) 



for z, z' < 0, whereas the latter two column vectors 
V\ and are obtained from and U^, respectively, 
by particle-hole conjugation (complex conjugation, inter- 
change of first and third, and of second and fourth rows, 
and the replacements £ — > — £ and ky — > — ky, see App. 
lAl). Here 



fcf2(£) = yfc| - fc| + e/hv^^, 



(13) 



up to corrections of order £^ , with v^z = fik^z (0) /m and 
''^iz = f^k\^z{^) /m. The coefficient p^{e) is the amplitude 



for normal reflection at the HS interface in the absence 
of the perturbations, see Eq. ([26]) below. Since there is 
no Andreev reflection in the absence of a magnetization 
gradient, one has |p(£)| — 1. The other coefficients do not 
have an interpretation in terms of scattering amplitudes. 
Detailed expressions for the coefficients pa and Tq- can be 
found in App. \K\ 

In the normal state (for Aq = 0), majority electrons in- 
cident onto the superconductor interface are transmitted 
with probability 



r^(k||) = 



^V^zVSt. 



(14) 



where vgz = {^/'m){kg — A:y)^/^. The use of the step- 
function model for the superconducting order parame- 
ter used in Eq. Q requires that the normal-state trans- 
parency of the HS interface is small for majority electrons 
as well as minority electrons .'SS In our model, this condi- 
tion is met if 



max(|ti;|, ug) ^ max(w^, w^). 



(15) 



and we will assume that this inequality is met in the 
results we present below. In our final results, we will 
eliminate w and vg in favor of the normal-state reflection 
amplitude r and the transmission coefficient T-^(k) given 
above. Expressions for the case of an ideal interface with 
transparency T-|-(k||) = 1, but still using the step- function 
model for A(r) at the interface, are given in App. [C| 



B. Perturbation: Changing magnetization and 
impurities 

The perturbation V describes the combined effect of a 
magnetization direction that slowly changes as a function 
of the coordinate z, as well as impurity scattering. We 
take the magnetization direction to be^ - 



m = (ei cos ip + 62 sin (p) sin 6{r) + cos 0{r), 
where the polar angle 



9iz) 



z/ld if z < 0, 
ff z > 0, 



(16) 



(17) 



while the azimuthal angle ip remains constant. Trans- 
forming to a spin coordinate system in which m points 
along 63 everywhere in space and expanding up to first 
order in the gradient dO/dz, one finds that the single- 
particle Hamiltonian H takes the form Hq described in 
the previous subsection, with the additional terrrPfil 



ih 
2m 

ih 
mid 



((72 COS If — (Ti sin If) 
(ct2 COS ip — ai sin ip) 



dO d d d6 
dz dz dz dz 
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The impurity potential is taken to be a sum of the form and 



(19) ^-ekiihW = 



where r^- is the position of the jth impurity and Uj its 
strength. The impurity strength Uj is related to the scat- 
tering cross section aj for majority electrons, 



"? 9 

m Uj 



(20) 



In the main text, we will consider the case that the impu- 
rities are randomly positioned with density nimp in the 
vicinity of the interface. Here it is important to point 
out, that riimp is the impurity density at the interface, 
and that there may be a different impurity density in 
the bulk of the half metal. (Our results below show that 
only impurities within a distance ~ from the HS 

interface contribute to Andreev reflection. Hence, only 
the impurity density near the interface enters in our fi- 
nal expressions.) The situation that all impurities are 
located precisely at the interface is discussed separately 
in appendix [B] 

Combining these two perturbations, we thus find 



V 



V^ + Vi 

-VI - V* 



(21) 



C. Scattering matrix 



There are two linearly independent solutions of the 
Bogoliubov-de Gennes equation ([I]) for each wavevector 
k|| . At a large distance from the HS interface (z ^ for 
the coordinate system used here), they can be taken to 
be of the standard form 



*ek„c(r) 



pik||-r||+ife^^(e)2 



g*kJ|T||-ifc^^(e)2; 



g''kI|T||-ifctz(-e)z 



(22) 



/ree(£;kj|,k||) 




/ 

I \ 



rho(£;kj|,k 

V ) 



r||-ifct^(e)2 





U/ 



(23) 



^ik;|T||-ifc^^(£)z 



/reh(e;kj|,k||) \ 










fhh(e;k'|,k|| 




which represents a state with specified incoming electron- 
like or hole-like quasiparticle boundary conditions. To- 
gether, the amplitudes Tee, ^ch, "The, and rhh define the 
scattering matrix S of the HS interface, 



5(k;|,kj[;£) 



rco(kj|,k||;e) rch(k'| , k|| ; e) 
nic(kn,k||;£) rhh(kn,k||;e) 



(24) 



The scattering matrix is unitary. Particle-hole symmetry 
gives the further constrainlj^ 



roo(k|j,k||,£) = rhh(-k'|, -k||, -e)* 
roh(kn,k||,e) = rhc(-kn , -k|| , -e)* 



(25) 



In the absence of the perturbation V, 5 is found from a 
solution of the Bogoliubov-de Gennes equation ([T]) with 
H = Ho, 



5(k;|,k||;e) 



k' kii 



Pt(k||;e) 
Pt(k||,-e)* 



(26) 



Inclusion of the term V leads to a shift S ^ S + SS of the 
scattering matrix, which can be calculated perturbatively 
by means of the Born series, 



frhc(kn,k||;£) = 

II II j/j 

n=0 



(e,k;|,h|V(gV)"|£,k||,e) 



(27) 

and similar expressions for Sr^c, Sr^h, and (Jrij^. Here the 
retarded scattering states |£,k||,e)^ and |e,k||,h)^ are 
solutions of the Bogoliubov-de Gennes equation ([I]) with 
Ji = Jio and particle-like or hole-like incoming boundary 
conditions, respectively. The advanced scattering states 
|e,k||,e)^ and |£, k||,h)^ are solutions of the same equa- 
tion, but with particle-like or hole-like outgoing bound- 
ary conditions. Explicitly, the wavefunction ^'^k,, e('') 
the electron-like scattering retarded state reads [compare 
Eqs. (p2l) and (pel)] 



*?k„e(r) 



o«k||-rii 



(28) 



V p;(-£)e«^^(--)^ J 
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whereas the wavefunction ^'e,k||,h(r) of the correspond- 
ing hole-hke scattering state is obtained by particle-hole 
conjugation, 



gik||-r|| 

( 

V 



(29) 
\ 



The wavefunctions of the advanced scattering states are 
(recall |p(k||,£)| = 1) 



ekii C 



(r) 



PT(-e)*ek„h(r). 



(30) 



Taken together, Eqs. (27)-(30l and the expressions ([10|- 
( 12 1 for the Green function Q contain all information rel- 



evant for a calculation of the Andreev reflection ampli- 
tudes to arbitrary order in the perturbation V. 



III. ANDREEV REFLECTION AMPLITUDES 

We now describe the calculation of the Andreev re- 
flection coefficients rhc(£; ky, ky) and rch(£; ky, ky) in the 
presence of a magnetization gradient and impurity scat- 
tering at the HS interface. All calculations are performed 
to first order in the magnetization gradient (i.e., to first 
order in Without impurity scattering, one then 

finds Andreev reflection amplitudes that are diagonal in 
the wavevectors ky and kj| and linearly proportional to 
the excitation energy e. This calculation was originally 
performed in Refs. I7I10I and the result is briefly sum- 
marized in Sec. IIII Al below. To first order in the im- 
purity potential, we find Andreev reflections that are off- 
diagonal in k|| and k'| . These amplitudes remain finite at 
e = 0, but vanish for k|| = k'| . In order to find the lead- 
ing contribution to the diagonal amplitudes rhe(£; k|| , k|| ) 
and reh(e;ky,ky) we thus need to go to second order in 
the impurity potential. The first-order and second-order 
calculations with respect to the impurity potential are 
given in Sees. IIIB and III C[ where we restrict ourselves 
to the relevant cases e = and e = 0, ky = kjj, respec- 
tively. 



A. No impurity potential 

To first order in the perturbation V, only the mag- 
netization gradient term T4n contributes to the Andreev 



reflection amplitude. One findJl^l 

i(re-*'^+*'^AoTf(k|| 



rho(e;kj|,k||) = (5k k; 



8fct. , ^t(k||) 



(31) 



in the limit of a low transmission r-^(k||) of the half- 
metal-superconductor interface, where we kept the sub- 
leading term proportional to r^(k||)^ because it appears 
with the (small) energy Aq in the denominator.!^ (The 
divergence of Eq. (31) for \e\ — > Aq is an artifact of the 
small- expansion; Equation (31) and Eq. (42) below 
are not valid for |e| in the immediate vicinity of Aq.) 

The essential feature of the Andreev reflection ampli- 
tude (31) is that it vanishes at the Fermi level, e = 0. 



As explained in the introduction, this is a consequence of 
the very special symmetries of the Hamiltonian ([ij in the 
absence of impurity scattering. We now show by explicit 
calculation that impurity scattering, in combination with 
the magnetization gradient, gives rise to Andreev reflec- 
tion into the half metal at zero energy. 



B. First order in impurity potential 

We now address the effect of impurities on the Andreev 
reflection amplitudes. We first consider a single impurity 
and postpone the discussion of the effect of a finite but 
low concentration of impurities Sec. |IIID[ We choose our 
coordinates such, that the impurity is located at position 
Ti — (0, 0, z\) with z\ < 0. The corresponding impurity 
potential then reads 



(32) 



With this choice of coordinates, the normal and Andreev 
reflection amplitudes are functions of the magnitudes fcy 
and fcy of the initial-state and final-state wavevectors ky 
and k'l and the angle between these vectors only. In par- 
ticular, with this impurity potential 7r-rotation symmetry 
around the z axis is restored, 



5(e;kj|,k||) = ^(e; -kj, , -k,. 



(33) 



Evaluating the zero-energy Andreev reflection ampli- 
tude The to first orde r in V,„ and first order in the impurity 
potential Vi of Eq. ( 32 ) , we find a nonzero contribution 



to the Andreev refiection amplitude at zero energy, 
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rho(0;kn,kji) 



2mie-**+^'^[e<^^'fct-r^(kj|)(cos(A;t.Zi) - e«^^^0 - K;,Tt(k||)e«^^^' sin(fc^^Zi)] 



(k„okj|), (34) 



where, as before, we have given the result to leading order in the transparency T-^ of the interface. [Unlike in the case 
of Eq. (31) there are no terms of subleading order in that come with small energy denominators.] The Andreev 
reflection amplitude The follows from the particle-hole symmetry relations ( 25 ) . The above result vanishes if the 
impurity is located precisely at the superconductor interface. In that case, the leading contribution to the Andreev 
reflection amplitude is of higher order in the interface transparency. We have listed the corresponding expressions in 

App.m 

The Andreev reflection amplitude (34) is odd under exchange of the initial and final momenta k|| and k'|. In 
particular, rhc(0; ky , k|| ) vanishes for Tcii = kj|. This can be seen from the following simple argument: For the 
calculation of the diagonal elements rhc(0; k|| , k|| ) to first order in and Vi, there is no difference between an 



impurity potential of the form ( 32 ) and a potential that is invariant under translations parallel to the superconductor 
interface, M(r) = {WxWy)'^Uio(z^ — z\). For the latter potential, the general arguments of Ref. ITUl applv. from which 
it follows that rhc(0; k|| , k|| ) = 0. 

Note that only impurities within a distance ~ 1/^^ of the minority electron wavefunction decay length from the 
interface contribute to the Andreev reflection amplitude The- This is consistent with the picture that the Andreev 
reflection process involves an intermediate (evanescent) minority electron state, which is then converted into a majority 
state via a spin- flip process enabled by the combination of the impurity scattering and the magnetization gradient. 
Impurities at a larger distance from the superconductor interface do not contribute to the Andreev reflection amplitude. 
Their contribution to observable quantities, such as the conductance of an HS junction or the Josephson current in 
an SHS junction can be calculated using standard approaches, see, e.g., Ref. |331 



C. Second order in impurity potential 



A nonzero contribution to the diagonal amplitudes rch(0; k|| , k|| ) can be expected in seconrf-order perturbation theory 
in the impurity potential Vi. Technically, it is most convenient to calculate the diagonal amplitudes roh(0; k|| , k|| ) from 
the leading order [first order in V\) results for the off-diagonal amplitudes roh(0; k'| , ky) of Eq. (34) and the first-order- 
in- 14 off-diagonal normal reflection amplitudes using the relation 



rhe(0;k||,k||) = - 



E 



3(0;k||,kj|)rhc(0;k||,kj|; 
roc(0;kji,kj|) 



(35) 



which is obtained upon combining unitarity of the scattering matrix, the particle-hole symmetry relations (25), and 
the symmetry relations (33). For the off-diagonal normal reflection amplitude roc(0; kn , ky) we find 



ree(0;k:,,k||) = 



4mi sin(fc.|.2Zi) sin(A;^^Zi) 



(36) 



to leading (zeroth) order in the interface transparency. To leading order in the interface transparency, the diagonal 
normal refiection coefficient roo(0;k||,k||) = —1. Because rijc(0; k|| , k'| ) vanishes for |k||| = |k'||, the summation in Eq. 

Ik:, 



(35) is dominated by terms in which jkii 



k^. The expression for the diagonal second-order- in- 14 Andreev 



reflection amplitude that we find upon substituting Eqs. (34) and (36) into Eq. (35) is too lengthy to be reported 
in full, but the result takes a simple form in the limit <^ min(w^,us) ^ niax(]u)|, ws), corresponding to a large 
mismatch of Fermi velocities in H and S, 



rhc(0;k||,k|| 



Abh?W.xWyKlvlldir 



[6(1 - e'^^^O - 6K^Zi + 3(^^Zi)2 + [n^zif 



(37) 



Equation (37) is valid for arbitrary impurity locations ri. Only impurities within a distance ~ 1/k^ contribute 
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to the Andreev reflection amplitude of Eq. (37 1. Note 



that the first order amplitude ( 34 ) and the second-order 



amplitude ( 37 1 are both proportional to the same power 



of the interface transparency r-^(k||). 



D. Finite density of impurities 

We now consider a finite but low density riimp of impu- 
rities in the immediate vicinity of the surface, for which 
the potential strength Uj is itself a random variable with 
zero mean and with variance (Uj). In view of the ap- 
plications of the next section, we are interested in the 
ensemble averages (rhc(0; ky , ky)) and (|rhe(0; ky , k||)p) 
with respect to the disorder at the interface, to lowest 
(first) order in niaip- The single-impurity results derived 
above are sufficient for this calculation, since interference 
effects between Andreev reflection processes that are en- 
abled by scattering off different impurities can be ne- 
glected to this order in the impurity density riimp. (They 
give a contribution proportional to nf^ip-) 

Thus proceeding, we find that the average reflection 
amplitude (rhc(0; ky , k|| )) is nonzero only if ky = ky 
(because translation invariance along the interface is re- 
stored upon taking the ensemble average), 



(rhc(0;k||,k||)) 



mm 



3 ^ 



where (a) = m? {u^) / nfi^ is the mean scattering cross 
section of the impurities, see Eq. (20). The mean square 
Andreev reflection amplitude is 

((|rhe(0;k;,,ky)|2) = ni„p(a)r^(kj|)r^(ky) 



(39) 



Both results are for the limit <C min(u^,us) <C 
max(|w|, lis). 

It is an interesting question, what the effect of impu- 
rities is on the symmetry of the superconducting cor- 
relations in the half metal. In the quasiclassical Green 
function approach, the superconducting correlations in 
the half metal are described with the help of the anoma- 
lous Green function /(if2;k, r), Vl being the Matsubara 
frequency, which is related to the Andreev reflection am- 
plitudes aJ^ 



/(zr!;k,r) 



2g-2|0|(r-e,)/;it,t. 

roh(if^;k||,k||) 
X < -rhe(-»f^;kj|,ky 




fc^ < 0, f7 > 0, 
fc^ > 0, f7 < 0, 
otherwise, 

(40) 



where k is a wavevector with fc = fc^. Without impuri- 
ties at the interface, one thus finds that the spin-active 



interface gives rise to correlations of predominantly odd- 
frequency s-wave type, with an anom alous Green func- 
tion / (X for small frequencies .l^^llll 'VVe find from Eqs. 
( 38 1 and ( 40 ) that the presence of impurities at the in- 



terface does not change this fundamental symmetry of 
the induced superconducting order, but it does change 
the asymptotic low-frequency dependence of /, / being 
proportional to sign (Tt) for small 57. This enhancement 
of odd frequency s-wave superconducting correlations at 
small frequencies is a unique signature of the impurity- 
assisted spin-triplet proximity effect in half metals. 



IV. APPLICATIONS 



A. Conductance 

The conductance of a half- metal- superconductor inter- 
face is given by the expressiorP^I^Sl 



|rhe(e = ey;kl,,ky)|2, (41) 



Without impurities, rho is given by Eq. (31) if there is 



a magnetization gradient perpendicular to the interface, 
and G is proportional to at low voltage,'^'^^ 



(38) G{V) 



64(A2 - (ey)2) 



(42) 



where T|(0) is the normal-state transmission coefficient 
of the clean HS interface for perpendicular incidence 
(k|| = 0), and we took the limit <C min(u^,us) <C 
max(|w|, vs)- 

For low bias, an impurity at position ri = Zje^ in- 
creases the conductance by a finite amount, which is 
found by substituting Eq. (34) into Eq. (41). Because 



of space limitations, we here report the resulting expres- 
sion for the limit ^ min(w^, us) ^ max(|w|, vs) only. 



SG = 



2eVir2(0)fc« , 

X [6 - 6e«^"' - 6K^Zi + 3(K^Zi)2 + (k^z;)' 



(43) 



where Ci is the scattering cross section of the impurity, 
see Eq. ( 20 ) . For a finite but small concentration rtimp 



of impurities near the superconductor interface one may 
add these contributions to the conductance, giving a fi- 
nite conductance at zero bias, 



G(0) 



437e2ni„,p(a)I^,W^^T^(0)2fc8 

34020772/1^9/2 



(44) 
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FIG. 2: SHS junction of length L with two spin-active inter- 
faces. 



where (cr) is the mean scattering cross section of the im- 
purities. At finite bias, the conductance is the sum of 
Eqs. (|42l) and (llil). 



The impurity-assisted contribution to the conductance 



(44 1 dominates over the conductance (42) of a clean HS 



interface in the limit of low bias voltages, 



Iml 



< 



* imp \ 



(45) 



where = ~^K^t'i/2 is the excitation gap for minority 
quasiparticles in the half metal, see Eq. ([8|. This condi- 
tion is independent of the normal-state transmission co- 
efficient T|(0) and the magnetization gradient ^, since 
the interface conductance G is proportional to T^{^)'^r^ 
for a clean interface as well as for an interface with impu- 
rity scattering. (Here we took the first term in Eq. ( 42 ) as 
the basis for our comparison, which is the leading contri- 
bution to the conductance of a clean interface in the limit 
of a small interface transparency r-^(O).) The fraction on 
the left hand side of Eq. (451 is the total cross section 



of all impurities within a layer of width adjacent to 
the HS interface, per unit interface area. 



B. Josephson current 

Next, we consider a superconductor-half-metal- 
superconductor junction, in which the central half- 
metallic segment is clean, with the possible exception of 
the presence of impurities with density riimp near the two 
interfaces. We choose coordinates such that the two in- 
terfaces are a.t z — ~L and z = 0, see Fig. [2] We consider 
the case that the magnetization gradient is perpendicu- 
lar to the superconductor interface at both interfaces, 
with equal azimuthal angles and that the magnetiza- 
tion gradients are equal in magnitude and opposite in 



direction. We restrict our analysis to the so-called "long 
junction" limit El <C Aq, where 

El = ^ (46) 



is the Thouless energy of the junction, and we consider 
temperatures in the range ^ <^ Aq. In this 

temperature regime, the Josephson current is quadratic 
in the Andreev reflection amplitudes,!^ 



I = - 



ieksT d 



^-TTkBTL/hv-f^-TrkBTL/hv!f^ 

kii.k: 



X Re Toh (iTrfceT; k|| , k'| )f he {iTrkBT; ky , ky ) ,(47) 

where /cb is the Boltzmann constant, T temperature, (f> 
the phase difference between the two superconductors. 
The amplitudes roh and fhe describe Andreev reflection 
at the interfaces at z = and z = — L, respectively. 

Each Andreev reflection amplitude that appears in Eq. 
(471 is the sum of two contributions: A contribution 



for the clean half-metal-superconductor interface, given 
in Eq. (31 1, and a contribution from impurity-mediated 



Andreev reflection. Upon taking the average over the 
ensemble of impurities, the impurity-mediated contribu- 
tions two factors r^h and fhe can be averaged separately, 
as they refer to two different interfaces. After taking the 
ensemble average, translation symmetry along the inter- 
face is restored, and the ensemble-averaged amplitudes 
(reh(«7i'fcBT; ky, k||)) and (fhe(j7rfcBr; ky , ky)) are zero ex- 
cept if ky = kn, see Eq. (38 1. We then find, again in the 



limit <C min(i;^,ws) ^ maxduil, ws), that 
^ _ eTrW,WymO)^EL 



-e -'^■b^/'El gin, 



(48) 



'l5fcBTr^(0) 120fcBrfc^ 8k^(T) 



Ao 



TTK] 



Upon lowering the temperature, Eq. (48 1 remains valid as 



long as kQT > El, and the Josephson current saturates 
for temperatures k^T ~ E^. 

Comparing the Josephson current for a clean interface 
and the impurity-assisted contribution, we find that the 
impurity-assisted contribution dominates if 



,(f7) / k^ 



> 



max(fcBT, El) 

Iml 



(49) 



where we again took the current in the clean case in the 
limit of a low transparency of the HS interface. Again, 
this condition is independent of and 1^, because the 
dependence of the Josephson current on these parameters 
is the same (proportional to T^l'^'^) for clean and dirty 
interfaces. 



V. CONCLUSION 

We have verified through explicit calculation that im- 
purities in the vicinity of a half-metal-superconductor 
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(HS) junction with a magnetization gradient perpendic- 
ular to the interface (spin-active interface) quahtatively 
change the dependence of the Andreev reflection amph- 
tude The on the quasiparticle energy e and, hence, the 
bias dependence of the interface conductance G and the 
temperature dependence of the Josephson current / of a 
superconductor-half-metal-superconductor (SHS) junc- 
tion. Without impurities, one has rho oc £ at low energy.^ 
In contrast, impurities give rise to a nonzero value of rhc 
for e = 0. Interestingly, the zero-bias conductance of a 
HS junction and the low-temperature Josephson current 
of an SHS junction increase upon increasing the impurity 
concentration at the interface. 

Although it was to be expected, based on general sym- 
metry considerations, that impurity scattering leads to 
a nonzero Andreev reflection probability |rhcP, we found 
the remarkable result that a nonzero amplitude (rhc) re- 
mains after taking an average over impurity locations 
and potentials. It is this coherent impurity-assisted An- 
dreev reflection that is responsible for the increase in 
the low-temperature Josephson current and the super- 
conducting correlations at low frequency that we predict. 
The nonzero average persists if an average over the Fermi 
surface is performed. 

There is an important difference between the impuri- 
ties in the immediate vicinity of the interface we consid- 
ered here, and impurities at a larger distance into the half 
metal. Only impurities in the half metal that arc within 
a minority-electron wavefunction decay length from the 
interface may enhance Andreev reflection. Within a qua- 
siclassical approach, such impurities are considered as 



part of the interface, and they are an integral part of 
the boundary conditions that have to be applied at the 
half-metal-superconductor interface. Impurities that are 
located farther into the half metal lead to scattering of 
quasiparticles before and after Andreev reflection, but 
such impurities are not involved in the Andreev reflection 
process itself. Their effect can be treated with standard 
methods from quasiclassics or the scattering matrix ap- 
proach, which combine normal-state propagation inside 
the half metal with the Andreev-reflection boundary con- 
ditions at the half-metal-superconductor interface.'^^Iini 
In any case, since the induced superconducting order is 
of s-wave type for both clean and disordered spin-active 
interfaces, impurities away from the interface will only 
have a minor effect on the superconducting order induced 
in the half metal. 
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Appendix A: Scattering states 

The full expression for the column vectors and 
appearing in the Green function is 



Vr 



Vi = 



\ 
/ 

/ ^^(^^y^-^k^Ae)z+K,{-e)z' \ 






The detailed expressions for the coefficients pa and t„ of Eq. (|12|) are 



-vsziviz - iv^z)£ + {vsz + inz + 2w)(iw^2 + 2wj)y/Al - £2 



Ti(e) 



VSz{viz + 


iv-i^)e - 


i4z + + 2w){-tv^, + 2w))^/Al - £2 ' 






2ie-'^vszVt,Ao 


VSz{viz 


+ iv-t^)^ 


+ (^^s. + i^z + 2w){-iv^^ + 2w))^Al - £2 


vszivu + 


iv^^)e ~ 


i^iz - ~ 2w)(iw^, + 2w))^Al - £2 


VSziviz - 


iv^^)e + 


^^Iz + i^i^ + 2w)(jwt^ + 2w)) V^o - ' 









(Al) 
(A2) 

(A3) 
(A4) 
(A5) 
(A6) 
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Appendix B: Impurities precisely at the HS interface 

If all impurities are located precisely at the HS inter- 
face {i.e., Zi = for all impurities), the impurity-assisted 
Andreev reflection amplitude vanishes to the order in the 
interface transparency that was required to derive Eqs. 



(34 1 and (37) of the main text. In this appendix we col- 



lect the main results of this article for the case Zi — 0. 

For impurities located at the HS interface it is not pos- 
sible to eliminate the parameters w and vqz in favor of 
the normal-state transmission T^(k||) of the HS interface 
only. In addition, we need the imaginary part Imr^(kj[) 



of the reflection amplitude for majority electrons incident 
on the HS interface from the half metal. 



r^(k||) 



~2iw - vsz + v^z 
2iw + vsz + v-fz 



(Bl) 



The transmission probability r(kj|) for majority elec- 
trons incident at the half-metal-superconductor inter- 
face with the superconductor in the normal state is 
r^(k) = l-|r^(k||)|2. 

For the first-ordcr-in-Vi contribution to the Andreev 
reflection amplitude one then finds 



mie"*'^+'''^K' jr^(k;)Imr^(k||) -f T^(kj|)Imr^(k;)] . s 

to leading order in the interface transparency. The normal reflection amplitude is, again to first order in the impurity 
potential and to leading order in the interface transparency, 

mi[rt(k|| )T^(K ) - 4Inirt(k|| )lmrUk[ )] 
rec(0;kj|,k||) = ^ -J^ (B3) 

For the second-order-in- diagonal Andreev reflection amplitude one finds 

wfe-*('^-'^)fc3(Imrt(k||))rt(k||)[Tt(k||)2 - 4(Imrt(ki|))2](5fc2 _ 3fc2) 

'■he(0;k||,k||) = ^ jy n2 2, .3 " ^ m 

where we took the limit <^ min(?;^,ws) <^ max(|ti;|, wg). In the same limit one finds that a single impurity at the 
HS interface increases the zero-bias conductance of a HS junction by the amount 

_ 8e2(a)T^(0)2[Imr^(0)]2fc| 

- 5257r^hKlll ■ ^^^^ 

Finally, for impurities located precisely at the superconductor interface with surface density ni„ip_2d, we find the 
Josephson current 



IbkBTTfiO) 120kBTk^ 2fc3(cr)nimp,2d(Imrt(0))[r^(0)2 - 4(Imr^(0))2 



Appendix C: Ideal interface 



(B6) 



Here, we analyze the case of an ideal half-metal-superconductor interface with perfect transmission, keeping the 
(non-self-consistent) step-function model for the order parameter A(r). Perfect transparency r-^(k||) = 1 at the 



NS interface is achieved by setting w = and vs = v^. In analogy to Eq. (34), the leading contribution to the 
impurity-assisted Andreev reflection amplitude, 

rhc(0;k'|,k||) = ' , j -2e''4--"^'fcj^fc^^[fc^^ cos(/c^^Zi) -f sin(fc^^Zi)] 

IdhW^Wyik^^ + nfjik^^ + iKiz){k\^ - ik[,)^v\^v^, 

+ 2e^^^^^k^znUk\, sin{k\,z,) - cos{k\,z,)] - e^^^^'+^^^'fc^.l^^^ - Sk^)} - (k,, kj,)] , (CI) 
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vanishes for ku = kn. For the first-order- in- normal reflection amplitude, we find in analogy to Eq. (36), that 



ree(0;k:i,k||) = 



-4mie'"^~*'^ k^zkLe'^^'^^''^'^'^''^ - (k^z cos(xfc^z) - k^z sin(xfc^z))(K^^ cos(a;fcL) - fcL sin(a;fcL)) 



(C2) 

The remaining results in this appendix are for the limit vg — corresponding to a large mismatch of the 

(minority) Fermi velocities in the half metal and the superconductor. For the diagonal contribution to rhc, we then 
find 



h'^WxWyldk^z TTVf 



(C3) 



Taking the average in impurity position in the half- metal ( 38 ) leads us to 



(rhe(0;k||,k||)) = 



(C4) 



As described in subsection |IV A[ we calculate the leading correction to the subgap conductance caused by a single 
impurity in the half metal. 



5G 



h TTll 



and find, after taking the average in impurity position, 

G(0) 



2e2 n{a)ni,r.pW^Wy k^ 



(C5) 



(C6) 



Finally, for the impurity-assisted Josephson current in the long junction limit for high temperatures <C ksT ^ Ao 
we find 



/ =- 



16ev^ W^Wyk^ ( 1 Y k 



L 



yn,^ 

V ki-l 



sin((/)) • e' 



27rfcBT 3nimp(a-) 



Ao 



(C7) 
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